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Abstract. Posing Kepler's problem of motion around a fixed "sun" requires 
the geometric mechanician to choose a metric and a Laplacian. The metric 
provides the kinetic energy. The fundamental solution to the Laplacian (with 
delta source at the "sun" ) provides the potential energy. Posing Kepler's three 
laws (with input from Galileo) requires symmetry conditions. The metric 
space must be homogeneous, isotropic, and admit dilations. Any Riemannian 
manifold enjoying these three symmetry properties is Euclidean. So if we 
want a semblance of Kepler's three laws to hold but also want to leave the 
Euclidean realm, we are forced out of the realm of Riemannian geometries. 
The Heisenberg group (a subRiemannian geometry) and lattices provide the 
simplest examples of metric spaces enjoying a semblance of all three of the 
Kcplcrian symmetries. We report success in posing, and solving, the Kepler 
problem on the Heisenberg group. We report failures in posing the Kepler 
problem on the rank two lattice and partial success in solving the problem on 
the integers. We pose a number of questions. 



1. Introduction 

Newton formulated and solved what we call today "Kepler's problem" - the prob- 
lem whose negative energy solutions are Keplerian ellipses. The essential backdrop 
to the problem is Euclidean three-space and its group of isometries and scalings. 
Can we pose Kepler's problem on an arbitrary metric space? What properties must 
the space have if Kepler's three laws, or ghosts of these laws, are to hold? 

In 'Foundations of Mechanics' ([!]), Abraham and Marsden formulate classical 
mechanics as dynamical systems on the tangent bundle or cotangent bundle of a 
Riemannian manifold, which they call "natural mechanical systems" . The Rie- 
mannian metric defines the kinetic energy. One must choose a potential energy. 
In order to formulate Kepler's problem on our manifold, we take this potential to 
be the fundamental solution to the Laplacian. We choose a point on the manifold 
to be our "sun," which is the delta function source of the fundamental solution. 
Following Galileo we assume that the choice of location of the sun does not matter: 
that is, we will assume that our space is homogeneous. 

The mildest Riemannian departures from Euclidean space are the spaces of con- 
stant curvature: the sphere and hyperbolic space. About a century and a half 
before 'Foundations,' Lobachevksy ([H]), one of the founders of hyperbolic geom- 
etry, posed the Kepler problem as a "natural mechanical system" on hyperbolic 
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space. Later, Serret ([E]) posed the Kepler problem on the sphere. We are grate- 
ful to F. Diacu for these references and his survey [Bj. Kepler's 1st and 2nd laws 
hold in each of the three constant curvature geometries: hyperbolic space, the 3- 
sphere, and the original flat Euclidean space. But Kepler's 3rd law fails for these 
non-flat geometries for the simple reason that they admit no continuous scaling 
symmetries, or "dilations." 

We argue that in order to even formulate Kepler's third law our metric space 
must admit dilations. But if a space admits dilations, and is not Euclidean, then 
it cannot even be Riemannian! (We sketch the proof of this fact below.) The non- 
Euclidean spaces which admit dilations are subRicmannian: they are the Carnot 
groups. The simplest Carnot group is the Heisenberg group. 

This observation brings us to our main problem: pose and solve the the Kepler 
problem on the Heisenberg group. We will pose it. We will not fully solve it. 
We will show that all periodic solutions to the Kepler problem on the Heisenberg 
group must lie on the zero energy surface, and that the problem is integrable when 
restricted to this zero energy surface. Such solutions are described in Figures [2] and 
[4] A modified version of Kepler's third law holds for the periodic solutions. 

To write down the Kepler-Heisenberg problem we must have an explicit expres- 
sion for the potential: it is the fundamental solution for the su&Laplacian on the 
Heisenberg group. Luckily, Folland ([5]) found such an expression in 1970. 

We will also attempt to pose and solve Kepler's problem on some lattices. Lat- 
tices almost admit dilations: we can scale by positive integer scaling factors, but we 
cannot invert these scaling factors. The integers form the simplest lattice. We will 
pose and solve a Kepler problem on the integers. Our 'solutions' are of a high school 
nature. (We apologize in advance if our treatment here embarrasses readers with 
any skill in numerical methods and discretization.) These solutions are indicated 
in Figure [5j 

We then try to pose and solve the Kepler problem on the integer lattice in the 
plane where we run into fundamental problems which lead us to believe that the 
very definition of a discrete dynamical system is not yet well formulated. The heart 
of this problem is that the differences of values of the lattice potential - that being 
the fundamental solution of the lattice Laplacian - are irrational. 

Dedication and Acknowledgements. This article is dedicated to the mem- 
ory of Jerry and in thanks for all his inspiration. We would also like to thank 
the GMC group, in particular, David Martin de Diego, Juan Carlos Marrero, and 
Edith Padron for inviting us to the summer school in 2011 outside of Madrid. The 
formulation of the Kepler-Heisenberg problem was inspired by talking with many 
of the participants at that summer school. 

2. Kepler's 3 laws in a metric space. 

Let's recall Kepler's three laws for the motion of planets around the sun. 
Kl. Planets travel on ellipses with one focus the sun. 

K2. Equal areas are swept out in equal times. This law is equivalent to the 
conservation of the planet's angular momentum about the sun. 

K3. Period-Size: The period T of an orbit and its size a (semi-major axis) are 
related by a universal monomial relation a? = CT . 

The Keplerian planet moves in a Euclidean space. Do Kepler's laws even make 
sense on a general metric space? If not, what restrictions must we impose on the 



3 



metric space in order to make sense of a particular law? We discuss what is required 
of our metric space in order to formulate the corresponding law. 

Kl. We can define an "ellipsoid" for any metric space X. Fix two foci S,F e X 
and a positive number 2a. Consider the locus of points x e X for which d(S, x) + 
d(x,F) — 2a. If this locus is to be a curve then the metric space must be two- 
dimensional, for example, a smooth surface. Kl requires then that X is a two- 
dimensional, or that its Keplcrian dynamics can be reduced to two-dimensions. 
Kepler's problem has been posed and solved satisfactorily on the two-sphere and 
on the hyperbolic plane as described in the introduction. Its solutions satisfy Kl. 

K2 is equivalent to conservation of angular momentum. Angular momentum is 
conserved if the kinetic and potential energies in Newton's formulation of the Kepler 
problem are invariant under rotations about the sun. This requires isotropicness: 
all directions in the metric space are the same, at least through the sun. The 
two-sphere and the hyperbolic plane enjoy rotational symmetry and hence K2. 

K3 is a scaling law. It is an immediate consequence of the fact that the New- 
tonian potential V(x) = c/\x\ is homogeneous of degree —1. This homogeneity 
implies the space-time symmetry x,t i— > Xx,X 3 ^ 2 t, which is to say: if x(t) solves 
Kepler's equation then so does \x(\~ 3 / 2 t). From one periodic solution x(t) we gen- 
erate a one-parameter family x\. The energy-period relation in K3 for this family 
follows from the scaling symmetry 

2.0.1. Kepler's third law for homogeneous potentials in Euclidean space. Any ho- 
mogeneous potential V on a Euclidean space enjoys a version of K3. Homo- 
geneity is a scaling symmetry: x n- Xx => V(x) n- V(Xx) = X~ a V(x). 
We try to extend the symmetry to time and velocities by a power law ansatz: 
(x, t, v) n- (Xx, X^t, X~ v v). Balancing the resulting scalings of potential and kinetic 
energies implies that v = a/2. The requirement v — dx/dt yields (3 = 1 + (a/2). 
We are led to the extended scalings 



In terms of curves j(t), which are maps from t to x-space, the scaling operation is 



One verifies that if 7 satisfies Newton's equation 7 = —\7V(-f) then so does 7a. 
(Use that W is homogeneous of degree —a — 1.) The scaling symmetry thus takes 
solutions of energy H to solutions of energy X~ a H. Now if 7 is periodic of period 
T and with typical size a, then 7^ is periodic of period X l) T = A 1+Q / 2 T and typical 
size Xa. We thus arrive at our modified K3: T 2 = Ca 2+a . 

2.0.2. Dilations. To have an analogue of K3, our metric space must, like Euclidean 
space, admit dilations. 

Definition 1. A dilation of the metric space (X, d) with center S € X and dilation 
factor X is a map 8\ : X — > X which fixes S and satisfies d(5\x, S\y) = Xd(x, y) for 
all x,y G X . We say that the metric space X admits dilations if there is a dilation 
of X with dilation factor X for each X > 0. 



(x,t,v) 1 y (Xx,X 1+a ' 2 t,X- a ' 2 v). 



7 (t) ^ lx(t) - A 7 (A-( 1+ «/ 2 >i). 



Spherical and hyperbolic metrics admit no dilations. K3 fails for both. 
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3. Keplerian SYMMETRIES. 
We will restrict ourselves to metric spaces X which 

• are homogeneous (1) 

• are isotropic (2) 

• admit dilations. (3) 

We will call these three properties the Keplerian symmetry assumptions. 

Historical Motivation. Newton's biggest victory was probably his derivation of 
Kepler's laws K1,K2, K3, from more basic laws: Galilean invariance, his equation 
F = ma, and the specific choice of force F as '1/r 2 .' From these laws he derived 
what we today call Kepler's differential equation q = — q/\q\ 3 and thence Kl-3. A 
subset of the Galilean group is the group of spatial isometries and this relates to 
homogeneity and isotropicness. Dilations, as discussed above, are included so as to 
get a version of Kepler's third law. 

We recall the formal definition homogeneity and isotropicness. Let Isom(X) 
denotes the group of isometries of X. Homogeneity asserts that Isom(X) acts 
transitively on X. Isotropicness asserts that Isom(X) acts transitively on the 
space of directions through any point Sel. The sphere and the hyperbolic plane 
are homogeneous and isotropic, but they do not admit dilations. 

Proposition 1. If a Riemannian manifold is homogeneous and admits dilations 
then it is a Euclidean space. 

Proof [sketch]. See Gromov [9], prop. 3.15. Gromov defines the metric tangent 
cone T p X of any metric space at any point p € X as the pointed limit (X, Xd) 
as A — > oo. This limit need not always exist, but it does exist for Riemannian 
manifolds and equals the usual tangent plane, with its induced Euclidean metric. 
If the metric d admits a dilation with scale factor A then (X, d) is isometric to 
(X, Xd). Letting A — > oo we see that such an X is isometric to its metric tangent 
cone T p X for all p. QED 

Consequently, if we insist on satisfying all three Keplerian symmetries (l)-(3) 
while also leaving the realm of Euclidean spaces, we must also leave the world of 
Riemannian manifolds! The simplest non-Euclidean metric space satisfying (l)-(3) 
is the Heisenberg group with its subRiemannian metric. 

4. Kepler's problem and the Laplacian 

Before formulating the Kepler-Heisenberg problem, we look into how the stan- 
dard Kepler problem fits within the framework of "natural mechanical systems" 
and thus how it generalizes to general Riemannian manifolds. This background 
will yield a straightforward way to place the Kepler problem in the Heisenberg 
context. 

The Hamiltonian for the standard Kepler problem on K 3 is 

h = ^(pI+pI+pD- j, 

where r = \J x 1 + y 2 + z 2 and a > 0. Why the 1/r potential? Perhaps the best 
answer is that U = is the fundamental solution for the Laplacian A on (Eu- 
clidean!) M 3 , i.e. the solution to AU — —Sq. (See [3] and references therein.) The 
choice of sign convention is due to the positivity of the operator — A. 
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The kinetic term in H is the principal symbol of the Laplacian, so we can write 

H(q,p) = ±a A -aA- 1 (*) 

where ox denotes the principal symbol of L, and A -1 (q) = K(q, 0), where K (x, y) — 
U(x — y) is the Green's function for the Laplacian (and where a = jr). This 
reformulation suggests that we can pose 'Kepler problem' as a Hamiltonian system 
on any 'space' A with a 'Laplacian' A. 

This prescription (*) for H leaves us with a number of puzzles. 

Problems. What is the cotangent bundle of an arbitrary 'space' A? Assuming 
we make sense of H as a function on the cotangent bundle of X, then what are 
Hamilton's equations on T* XI Can we ever compute the fundamental solution 
A^ 1 of our Laplacian? 

All these questions have answers in the Riemannian case. The principal symbol 
has the coordinate expression 

<r A (p) = (q)piPj 

— it is the standard cometric of kinetic energy. The fundamental solution of the 
Laplacian has been explicitly computed for hyperbolic n-space, so we have a hy- 
perbolic Kepler problem. 

If A is a compact manifold without boundary, then the fundamental solution 
A^ 1 docs not exist for topological reasons. For example, we cannot have a single 
gravitational source on the sphere. There must be an opposing sink elsewhere on 
the sphere. To formulate the Kepler problem on the sphere, one places the sink 
antipodally to the source. See [5] or [TS] for a precise formulation. 



5. Kepler's Problem on the Heisenberg Group! 

5.1. Heisenberg Geometry. Consider M 3 with standard x,y,z coordinates, en- 
dowed with the two vector fields 

x _ d 1 d _ d 1 d_ 

dx 2^ dz dy 2 dz 

Then X, Y span the canonical contact distribution D on M 3 with induced Lebesgue 
volume form. Curves are called horizontal if they are tangent to D. Declaring 
A, Y orthonormal defines the standard subRiemannian structure on the Heisenberg 
group and yields the Carnot-Caratheodory metric dsjjj = dx 2 + dy 2 . Geodesies 
are qualitatively helices: lifts of circles and lines in the xy-plane. The horizontal 
constraint implies that the 2-coordinate of a curve grows like the area traced out 
by the projection of the curve to the xy-plane. See Figure [I] and Chapter 1 of [H] . 
The Heisenberg (sub)Laplacian is 

A = A 2 + Y 2 , 

a second order subelliptic operator, and the only correct choice for 'Laplacian' on 
the Heisenberg group. We have 

d 

{X,Y] = -=:Z 
dz 

and [A, Z] = [Y, Z] — 0. There are the commutation relations of the Heisenberg 
Lie algebra, hence the name. The Heisenberg group H is the simply connected Lie 
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Figure 1. A Heisenberg geodesic 

group with Lie algebra the Heisenberg algebra and is diffeomorphic to K 3 . In x, y, z 
coordinates the Heisenberg group law reads 

{xuVi,Zi) ■ (x 2 ,V2, z 2 ) = (xi + 2:2,2/1 + 2/2, zi + z 2 + \(x-iy 2 - x 2 yi)). 
Left multiplication is an isometry and the vector fields X, Y are left invariant. 

5.2. The Heisenberg Kepler Problem. Folland ([5]) has derived an explicit 
formula for the fundamental solution for the Heisenberg Laplacian! It is 

Here a = 2/ir. Let p Xl p Vl p z be the dual momenta to x,y,z so that together 
x,y, z,p x ,py,p z form canonical coordinates on T*H. Then 

Px = Px ~ \vPz, Py = Py + \xp z 

are dual momenta to X, Y, and 

K=\{P 2 x +P Y ) = \a A 

is the Heisenberg kinetic energy, given canonically by the cometric. (See Chapter 
1 of [Mj.) K generates the subRiemannian geodesic flow on the Heisenberg group. 
We see that Keplerian dynamics on the Heisenberg group are the Hamiltonian 
dynamics for the canonical Hamiltonian 

H = K — U. 

There is no explicit formula for the Heisenberg subRiemannian distance function 
\\(x, y, z)||h := d sr ((x, y, z), (0, 0, 0)), measuring the distance from a point to the 
origin . So the mix of K and U - of geodesic and subLaplacian - is quite interesting 
and it is rather remarkable that we can write down the Hamiltonian in closed form. 

The dilation on the Heisenberg group is 

h(x,y, z) = (\x,Xy,X 2 z). 



7 



Like the subRiemannian distance, the function p is positive homogeneous of degree 
1 with respect to this dilation. Since the Heisenberg sphere is homeomorphic to 
the Euclidean sphere, the standard argument which shows that any two norms on 
W 1 are Lipshitz equivalent shows that p and || • ||h are Lipshitz equivalent: there 
exist positive constants c, C such that cp(x,y 1 z) < \\(x,y, z)\\a < Cp(x,y,z) for 
(x,y,z)^0. 

Following the procedure described in Section |2.0.1[ we find that if a curve 7 
solves Newton's equation 7(i) = VJ7(7(i)), where V denotes the subRiemannian 
gradient, then so does 

7a(<) :=<5a(7(A- 2 <)). 



Then given a periodic orbit 7 with period T (see Section 5.5|, we get a family of 
periodic orbits 7a with periods X 2 T. Choosing a suitable notion of the 'size' a of a 
periodic orbit yields the Heisenberg version of Kepler's third law: 

T 2 = Ga\ 

The isometry group of the Heisenberg group is generated by translations and 
rotations. The translations denote the action of the Heisenberg group on itself by 
left multiplication. These project to translations of the xy-plane. The rotations 
form the circle group of rotations about about the z axis. In addition we have 
the discrete 'reflection' (x,y,z) <— > (x,— y, — z). Translations act transitively: the 
Heisenberg group is homogeneous. Rotations act transitively on (allowable) direc- 
tions: the Heisenberg group is isotropic. Thus the Heisenberg group enjoys the 
three Kcplerian symmetry properties. 

5.3. Hamiltonian Dynamics. The dilation on phase space T*H is 

<5 A : (x, y, z,p x ,p yi p z ) H> (Xx, Ay, X 2 z, X' 1 p x , X^p y , X~ 2 p z ). 

This is generated by the function J = xp x + yp y + 2zp z , which satisfies J = 2H. 
When H — 0, J is a first integral. Note that d\: H n- X~ 2 H. 

Now change to cylindrical coordinates (r,8,z) on H. We have the induced con- 
jugate momenta p r = {xp x + yp y )/r and pg — xp y — yp x . Our Hamiltonian is 



1 2 , lfPo , 1 \ 2 -2 

2Pr+2\— +2 r Pz) ~ a P ■ 



H=2 

Note that this does not depend on 9 due to rotational symmetry, and the corre- 
sponding angular momentum pg is conserved. 

On the smooth submanifold of phase space {H — 0}, we have three (independent) 
conserved quantities H,pg, and J, and a theorem of Arnold (see [2]) says that our 
system is intcgrable by quadratures here. See Figure [2] for approximations of orbits 
which exhibit this integrable behavior as well as the helical Heisenberg geometry. 
For this reason, we will mostly focus on the H = case. This is especially justified 
in light of the following. 

Lemma 1. Periodic orbits must have zero energy. 

Proof. If j(t) = (x(t),y(t),z(t),p x (t),p y (t),p z (t)) satisfies 7(0) = j(T) for some 
t = T, then J = xp x + yp y + 2zp z is also periodic. But we know the time derivative 
of J is constant, given by J = 2H. Since J cannot be monotonically increasing nor 
decreasing, we must have J = 2H = 0, so H = 0. QED 
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Periodic orbits exist and the existence proof forms part of C.S.'s thesis - see 
Section 5.5 below. We will momentarily report progress with integration of the 
H = system, but first we gather other dynamical results. 

Proposition 2. If H < then any solution is bounded. 

Proof. Suppose H = —h where h is positive. Then K — U — —ft, so 

U = K + h> h, 

since K is always non-negative. Then a solution (x(t),y(t), z(t)) in configuration 
space must satisfy 

where a and h are positive constants. QED 
Proposition 3. The only solutions in the plane z — are lines through the origin. 
Proof. The equations for z and 9 satisfy the relation 

z = \r 2 6. 

For a path lying in the plane 2 = 0, this implies either r = or 6 = 0. In the first 
case, the path is trivial. In the second, it lies on a line through the origin. Such a 
curve may be parametrized by 

lit) = (di 1 /*, paf 1 / 3 , 0, id*" 1 / 2 , &t~ 1/2 > °)- 
It is easy to verify that the desired equations are satisfied, and that H — 0. QED 

Proposition 4. The only solutions constant in configuration space are 

7(t) = (0,0,M,0,-fjt). 

Proof. This is an easy calculation. Note that such solutions are unbounded in phase 
space, and satisfy H < 0. QED 

Next, we explicitly integrate the equations of motion on a codimension 3 sub- 
manifold, and recover conies reminiscent of the Euclidean Kepler problem. Consider 
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the smooth submanifold N = {z = p z = pg — 0}. This submanifold is invariant 
under the dynamics, since z = p z = pg = on N. The Hamiltonian is 

which has the form of a classical central force problem in the plane. Fix an energy 
level H\ N — h. Then since p r — r, we can explicitly solve for r(t) as follows. 

Proposition 5. On N, r(t) traces out a hyperbola if h > 0, an ellipse if h < 0, 
and a parabola if h = 0. 

Proof. The Hamiltonian may be rewritten as the simple ODE 
Assume temporarily that h ^ 0. Integrating, we find 

which may be rewritten r 2 — 2ht 2 = — f • Since a > 0, this curve in the t, r-plane is 
an ellipse for h < and a hyperbola for h > 0. 
If h = 0, we find that 

and thus r 2 = y/8a t. QED 
We conclude this section with the following conjecture: 

Conjecture. There is an open set of initial conditions whose orbits are asymptotic 
to helices. 

This behavior is suggested by numerical experiment and by the fact that U and 
its derivatives tend to zero as orbits tend towards oo. 

5.4. Integration of the case H = 0. We now focus on the H = case and reduce 
the integrability of the equations of motion to the parametrization of a family of 
degree 6 algebraic plane curves. 

Let H = jj. Then integral curves for H are the same as geodesies for the metric 
Uds^. When H = 0, this is the same as the metric (H + U)ds^, whose geodesies 
correspond to integral curves for H, according to the Jacobi-Maupertuis principle. 
Thus, the flow of H is the same as the flow of H up to reparametrization on the 
hypersurface {H = 0} = {H = 1}. 

A short calculation shows that both J and pg Poisson commute with H. (Recall 
{H, J} = 2H.) This demonstrates the scale invariance of H ; 6\ : H i-> H . More 
importantly, we have three independent quantities conserved by the flow of H. 
Thus, we have an integrable system on {H = 0} = {H = 1} . 

Change our third coordinate z v = z/r 2 . We have conjugate momenta 
p r ,pg,p v . In these coordinates, we have J = rp r (as in the Euclidean case) and 
H = H(pg, J,v,p v ). On the submanifold {H = 0}, we have H = 1. Also, the initial 
conditions determine the constants J and pg. Thus, given initial conditions, H is 
a function of v and p v only. We arrive at the following result. 
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Figure 3. Curves in the v,p„-plane corresponding to J = 3, pe — 
1 (left) and J = 0, p g = 1 (right). 

Proposition 6. When H = 1, any solution must project to an algebraic curve in 
the v,p v -plane. 

These curves are naturally degree 10 but can be reduced to degree 6 by changing 
variables. Examples are shown in Figure [3] If we can parametrize these curves, we 
should be able to bootstrap up to find explicit solutions. 

5.5. Periodic Orbits. Despite the fact that the H = case is integrable, we have 
not been able to explicitly solve the equations. However, we know that periodic 
orbits exist. 

Take L = K + U as our Lagrangian and impose the horizontal constraint z = 
\xy — \yx. Then any trajectory 7 must lie on the zero set of the function 

G = \xy — hyx — z. 

The calculus of variations tells us that if 7 : [0, T] — > H is a minima of the action 
functional J Q Ldt which also satisfies our constraint, then there exists a scalar 
A = A(t) such that 7 is a minima of the modified action functional 

A(7)= / L x {t,-y,j)dt, 
Jo 

where we have written L;\(t, 7, 7) = L(t, 7, 7) — A(i)G(f, 7, 7). Setting the first 
variation of A equal to zero and integrating by parts yields the Euler-Lagrange 
equations: 

x = -Ay - \\y - 2ax(x 2 + y 2 )p~ 6 
y = Xx+ \\x - 2ay(x 2 + y 2 )p~ 6 

When A = p z we find that these agree with Hamilton's equations. 

We are confident that the direct method in the calculus of variations applied to 
A(-f) will yield a proof of the existence of periodic orbits. One works in the Hilbert 
space 7? 1 (5' 1 ,IHI) and requires that admissible curves are horizontal and satisfy the 
symmetry conditions 

l(t + T/3) = R 27l/3l (t) (SI) 

and 

z(t + T/2) = -z(t), (S2) 



Figure 4. Projection of an orbit to the xy-plane (left) and z- 
coordinate over time (right) 



where 

2 2 u 

2 2 u 

1 

Any admissible curve is therefore necessarily periodic, with additional symmetry. 
A suggestive approximation of such a curve is shown in Figure [4j 

The idea is to choose a minimizing sequence j n of curves in this space, and show 
that they converge within the space to some 7*. Applying elementary analysis 
and the principle of symmetric criticality shows that 7* must minimize the action, 
thereby satisfying the Euler-Lagrange equations. A central difficulty lies in proving 
that 7* does not pass through the singularity at the origin. A full existence proof 
is expected to appear in the thesis of C.S. 



#2tt/3 — 



5.6. A failure of reduction. Newton reduced his two-body problem to the Kepler 
problem in Euclidean space. There is no analogous reduction for the two-body 
problem on the Heisenberg group, nor is there for the two-body problem on the 
sphere or in hyperbolic space. We discuss the geometric roots of this failure. 

We begin by writing down the Heisenberg two-body problem. Let qx,q% € H = 
K 3 denote the positions of two bodies moving in the Heisenberg group H. Let their 
masses be 7711,7712. Their individual kinetic energies are 

1 

2771,; 



where P l x ,Py are the horizontal momenta of each body, as in Section 5.2 The 
Heisenberg two-body problem is defined by the Hamiltonian 

H = K\ + K2 - Kmim2U{q^ 1 q2), 

where k is the Gravitational constant and U is Folland's fundamental solution. H 
is a Hamiltonian on the cotangent bundle of H x H, and is invariant under the 
(cotangent lift of the left) translation (qi,qi) H> (391,392), 
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We know of two derivations of Kepler's problem (on Euclidean space) from New- 
ton's two-body problem. We will call these the 'algebraic' and the 'group-theoretic' 
derivations. The 'algebraic derivation' begins with the equation F = ma for each 
body. Divide the equation for each body by its mass to get equation for the ac- 
celeration q\ of each body's position vector q^. Subtract one equation from the 
other to obtain the ODE of Kepler's problem, q = — aq/\q\ 3 , for the difference vec- 
tor q = qi — q 2 . The 'group theoretic derivation' depends on the conservation of 
the total linear momentum, the invariance of Newton's mechanics with respect to 
Galilean boosts, and the abelian nature of the translation group. If P is the total 
linear momentum and M the total mass, we boost by the velocity —P/M to get 
to a new representation of the same dynamics in which the total linear momen- 
tum is zero. Then we reduce by translation at the value by placing the center of 
mass at the origin. Finally, we compute that each mass separately satisfies Kepler's 
equation with the origin - the center of mass - now playing the role of "sun" . 

The algebraic derivation fails on the Heisenberg group because the 'difference 
vector' gi(t)~ 1 g2(t) of two Heisenberg geodesies is not a Heisenberg geodesic. Why 
is this lack of being a geodesic a problem? Set the Heisenberg Gravitational constant 
n = so the two-body problem reduces to two uncoupled Heisenberg geodesic 
problems. Play the algebraic game. Our 'difference vector' does not satisfy the 
Heisenberg geodesic equations or any other pretty Hamiltonian equation. But in 
the Newtonian-Euclidean case, the difference vector travels like a free particle, i.e., 
moves in a straight line - as it should with a = in Kepler's problem. Things will 
just get worse for n ^ 0. 

The failure of the group theoretic derivation goes a bit deeper and is perhaps 
more enlightening. What is a 'Galilean boost' for an arbitrary Lie group? We choose 
some 'translation velocity' £ and multiply elements xq by exp(t£) . Euclidean space 
enjoys the wonderful property that exp(t£)xo — xq + i£ describes free motion; it is 
a geodesic. This assertion is decidedly false for the Heisenberg group: the orbits of 
(left or right) translates of one parameter subgroups are not Heisenberg geodesies. 
As a result, applying a boost to a solution (qi(t), q 2 (t)) to the Heisenberg two-body 
problem will not yield a solution. There is a conserved total 'linear momentum': 
the momentum map for the (left) translation action. But we cannot use it to 
'Galilean boost' the 'center of mass velocity' down to zero. Even if this total linear 
momentum were initially zero, we still seem to be stuck. The non-Abelian nature 
of the group appears to block us from writing the reduced Hamiltonian at zero as 
a Kepler Hamiltonian on the 'diagonal group' of elements q = g-f q%. 

In spherical and hyperbolic geometry, reduction of the two-body problem to the 
Kepler problem fails for similar reasons. See [7]. In the spherical case, Shchepetilov 
[IB] used the Morales-Ramis theory to prove that the two-body problem in these 
two geometries is not meromorphically integrable. 

Question. Is the two-body problem on the Heisenberg group non-integrable? 



6. Kepler's Problem on a Lattice. 

Lattices admit one-sided dilations: we can scale a lattice L by a positive integer 
c and land back in the lattice, stretching all distances by c. They admit Laplacians. 
So we might be able to begin to investigate Kepler's 3rd law on L. 
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What are Newton's equations on L? Since we must hop from lattice site to 
lattice site, we must choose our time variable t to be discrete: 

t= ...,-1,0,1,2,.... 

A 'solution' to Newton's equations will then be a 'discrete curve' 

7 : Z — > L, L our lattice, 

satisfying a difference equation which mimics Newton's equations. In 1st order 
Hamiltonian form these equations should resemble 

dj 

W] = " 
|j~WM<». 

where the differential is the discrete difference operator 

|p7(* + l)-7(*), 

and where V : L — > R is our potential. The standard interpretation of VV is in 
terms of its differential 

dV(l) :Ei^R, £ eh, 

where E = E^ is the set of edges (chosen lattice generatos) leaving the lattice site 
£, and where 

dV(£)(e) = V{£') - V(£), e = [£,£'] an edge. 
Then we can rewrite our Newton difference equations as 

7(t + l)=7(*)+P(t) (4) 
P (t + l)=p(t)-dV( 1 (t)). (5) 

What is the momentum, p(t)l We add to it dV(j(t)), so it must lie in the same 
space as dV, which is 

E* = E^ = real valued functions on E^. 

This 'cotangent space at P is a vector space isomorphic to R d , where d is the degree 
of a vertex: the number of edges leaving £. Good. Now, how do we add p{t) e E* 
to a lattice site £ = j(t) e L in order to get a new lattice site *y(t + 1) as in the 
1st Newton equation? We seem to be missing the 'mass matrix' or 'cometric' of 
mechanics. 

Definition 2. A lattice cometric is a 'non-trivial 7 map 

M : E* -> L. 

With this tentative definition we can now try to write down 'Newton's equations' 

7(t+l)= 7 (t)+M(p(i)) (6) 
p{t + l)=p{t)-dV( 1 {t)), (7) 

which define a discrete dynamics on the phase space L x E* . The resulting dynamics 
have some vague relation to the corresponding formal Hamiltonian 

H{£,p)= l -pMp + V{£), 
but we aren't sure how to interpret the term pMp. 
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Figure 5. A Kepler orbit in integer phase space 

7. Kepler's problem on Z. 

The Laplacian A on Z is given by A/(n) = f(n + 1) — 2f(n) + f(n — 1). 

Exercise. Show that U(n) = — ^\n\ is a fundamental solution for the Laplacian 
on Z with source at S = 0. 

Take the Kepler constant a = 2 so that the 'Newtonian potential' V = —all is 
V = \n\. The Hamiltonian is 

H(n,p) = ^p 2 + \n\. 

Since \n\ = sgn(n) is the sign of n, 1 if n > and — 1 if n < 0, we hnd that 
with this choice of a the discrete gradient is integer valued. We get a good discrete 
dynamical system. Newton's equations (in 1st order form) become 

n{j + 1) = n(j) +p(j) 

P(j + !) =pU) - sgn(n). 
A solution is depicted in Figure [5] 

At each iteration, the 'lattice momentum' p decreases by one as long as n > 
and increases by 1 as long as n < 0. (We have to make a choice at n = 0; above we 
chose sgn(0) = —1.) Note that as long as the initial condition p(0) is an integer, it 
remains an integer, and we stay on the lattice! 

Because of this happy coincidence with p's evolution, we did not need to worry 
about where p lived. It is a real number that happens to evolve to stay integral. 
Our momentum space is not E* = M 2 . (There are two directions, right and left, 
on the lattice, hence the dimension 2.) We also did not need to choose a 'lattice 
cometric' M. If any choice was made, it seems to have been '1' as written in the 
Hamiltonian. The happy coincidence does not happen when we go up to the rank 
2 lattice. 

7.1. Kepler on the rank 2 lattice. The rank 2 lattice is Z 2 with elements written 
I = (n,rn). As a metric space, we use the distance 

G?((n, m), (n' , m!) = \n — n'\ + \m — m'\. 
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Figure 6. Values of the Green's function gonZ 2 , taken from [TU] 
The Laplacian is 

#: d(^')=l 

Let £/ denote the fundamental solution, this being the 'most bounded' solution to 
AU = So, where So is the lattice delta function corresponding to placing the sun 
at the origin. (There is a lattice Liouvillc theorem, so U can be made unique up 
to an additive constant.) There is no closed form expression for U. However, any 
particular value of U can be computed recursively. Indeed, the fundamental solution 
is a well studied object with applications to the theory of electrical circuits ([!]), 
solid state physics, and quantum mechanics. Some values of the lattice Green's 
function are reproduced in Figure [6] from [10]. (We thank the brothers Hollos for 
permission to reproduce their table.) 
We write the formal Hamiltonian 

H(e,p) = ~pMp-aU(£) 

and derive Hamilton's equations 

£(t+l) = e{t)+Mp(t) 
P (t + l)=p(t)+adU(£(t)) 
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where 

M : R 4 = E* -> Z 2 . 

The vector p is a 4- vector with components (p up , Pdown, Pright, Pie ft) correspond- 
ing to the 4 edges, which are the 4 directions of motion, through each vertex. We 
have, for example dU(£) up = U(£ + e 2 ) — U(£) where e-i = (0, 1) represents motion 
in the 'up' direction. The second Hamilton equation makes sense. 

When we try to parse the first Hamilton equation we get stuck. What do we 
take for M : R 4 -> Z 2 ? We require M to be non-constant. Certainly M will not be 
continuous! Ideally M is 'linear': 

M(kp) = kM(p), fceZ, 

but this is probably not possible in any reasonable sense. One possibility for 
M is to argue that there is a 'canonical' projection n : M 4 — > M 2 , for example, 
i Pdown i Pright-) Pie ft ) — 2(Pu P + Pdown, Pright + Pie ft), and a canonical embed- 
ding of our lattice as Z 2 cK 2 . Then choose M(p) to be the lattice point closest to 
n(p). This leaves us to worry about what to do if n(p) is midway between lattice 
points. Flip a coin? 

We are stuck and look forward to some of our readers unsticking us. 

7.1.1. Euler- Lagrange formulation. Wc can make a bit more sense of the Euler 
Lagrange version of lattice dynamics. Fix a positive integer T, the 'time of flight,' 
and initial and final vertices, vq,v\ G Z 2 . There will be two formulations. In both, 
we consider discrete paths 7 : {0, 1, . . . , T} — > Z 2 which join v to v\ in time T, and 
we minimize an 'action functional' A among all such discrete paths. 
Version 1: Minimize the action 

T 

^(7)=E4(l^7Wli) 2 + ^(7W)} 

among all discrete paths 7 joining v to V\ in discrete time T. 

Here I^j7(i)|i = d(j(t + 1), 7(*)), so hah" of its square represents kinetic energy. 

Version 2: Call a discrete path 'continuous' if either d(j(t + l),-y(t)) = 1 or 
j(t+l) = j(t). Minimize the same action as Version 1, but now over all continuous 
paths. (In this case the kinetic term ^{\-^j{t)\i) 2 is either 1/2 or at each time 
step.) 

We are guaranteed a solution to Version 2 exists since there are only a finite 
number of 'continuous' paths joining v to v\. We suspect that if we move too fast 
the kinetic energy becomes too large, so that Version 1 is 'coercive' and one can 
argue that again there are only a finite number of paths that matter. 

It seems doubtful that any decent Euler-Lagrange type difference equation 'dy- 
namics' will result from either principle. Indeed, take the case a = of a 'free 
particle' on the latttice, and take vq = (0,0), v\ — (n,m), n > m > 0. There are 
(n + l)m shortest paths from v to v\. Just draw box-paths, always moving either 
right or up. Their lengths are all n + m = d(vo,Vi). If T < d(vo,Vi) then there 
are no paths connecting the two points. If T = d(vo,v\) then their actions are all 
\T. If T > d(vo, vi) the action remains the same; we just stay still for the requisite 
times T — d(vo,vi). This means either (i) all points v\ — (n,m) with nm ^ are 
conjugate to Vq, or (ii) that there is no good 'free' dynamical equation, so likely no 
good Euler-Lagrange equations in general. 
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8. Quantum Mechanics to Classical Mechanics on Cayley Graphs? 

By a graph here we mean the usual combinatorial collection of vertices and edges. 
We write T for the set of vertices and view T as 'configuration space.' The graph 
Laplacian is the operator A : ^(r) — > ^(r) defined by 

Af(v) = ]T (/(</) - /(«)). 

v':[v,v'\ an edge 

If T is finite there will be no fundamental solution; that is, there is no solution to 
AUs — Ss where where Ss is the discrete 8 function centered at the sun: Ss(S) = 
1, 6s{v) = 0, v ^ S. If r is finite, a necessary condition for the solvability of AV = f 
is = 0, which will fail for / = Ss- 

Regardless of whether or not T has a Green's function, it has plenty of potentials, 
meaning functions V <G ^(r). Consequently for each choice of Planck's constant h 
we have a Schrodinger operator: 

h 2 A + V :£ 2 (T)^£ 2 (T). 

There is a large active field of graph Laplacians and quantum mechanics on graphs. 
There is undoubtedly a theory of quantum mechanics on T. 

Challenge. Don't you think this quantum mechanics ought to have a classical 
limit? If 'yes' then please answer: what are the correct Newton's equations for an 
arbitrary potential, on an arbitrary graph? 

8.0.2. Cayley graph of a group. Let T be a finitely generated group and ei, . . . , e<j <E 
T be a fixed set of generators for T (so every element of T is a product of the e, 's 
or their inverses). Form the graph whose vertices are the elements x € T and for 
which two vertices x,y <G T are joined by an edge if and only if either y — xei or 
x = yei for some generator ej. Count each edge as having length 1. Define the 
distance between points x and y in T to be the minimum of the lengths of the paths 
joining x to y. This distance is always an integer, since the length of a path is just 
the number of edges it contains. 

In this representation, the 'Lie algebra' of the Cayley graph will be the tangent 
space at the identity: the disjoint union of d copies of Z. Alternatively, it is the 
subset of 7L d consisting of vectors for which all but one component is zero. 

Example: Lattices. Take r = Z 2 to be the lattice of integers in the plane, with 
standard generators e\ — (1,0), e 2 = (0, 1). Then the Cayley graph of Z 2 realized 
as above has the vertices of a standard infinite sheet of graph paper in R 2 . Its 
Lie algebra consists of integer points on the cc-axis unioned with the collection of 
integer points on the y-axis. 

8.0.3. Kepler symmetries of Cayley graphs. Every Cayley graph satisfies Kcplerian 
symmetry property (1) of being homogeneous since F acts on itself on the right 
by isometries. View the generators as the 'directions.' Then if the automorphism 
group of the group T acts transitively on its generating set ei, . . . , e<j, the metric is 
isotropic; it satisfies Kcplerian symmetry property (2). Finally we can send ei to 
e\. In some instances this defines a group homomorphism of T into itself. Then 
the Cayley graph admits one-sided dilations and so satisfies (3). The examples we 
know of groups whose Cayley graphs satisfy (1), (2) and (3) are the lattices Z d , the 
lattices in nilpotent groups, and the free group on d generators. In the continuous 
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case, we know how to derive a Kepler's third law from the Keplerian symmetry (3). 
Is there an analogous construction in the discrete case? 

8.0.4. Full Disclosure - R.M.. I have little interest in any kind of graph for its own 
sake. I am not a combinatorist, nor a discrete group theorist! 

In contrast to the dozens of books that Jerry wrote in his life, I have mustered 
the courage and stamina to write a single book in this life. (Jerry continues to 
amaze.) In that book I devoted a chapter to trying to understand one of the big 
ideas of Gromov in his paper 'On Groups of Polynomial Growth . . . ' ([8 1) in which 
he used subRiemannian ideas to solve a problem in discrete group theory. Consider 
a discrete finitely generated group V . Select some generators and form the group's 
Cayley graph. We say the group is 'of polynomial growth' if the number of vertices 
of the Cayley graph lying inside a ball of radius R is bounded by a polynomial in R 
as R — > oo. (If r is of polynomial growth with respect to one set of generators, it 
is of polynomial growth with respect to any other set of generators.) The lattices, 
and the integer lattice in the Heisenberg group are examples of groups of polyno- 
mial growth. More generally, the lattices in any Carnot group are of polynomial 
growth. The free group on 2 generators in not of polynomial growth: its balls 
have exponential growth, roughly 3 R . There is a notion of a group being 'virtually 
nilpotent,' and it was known that virtually nilpotent implies polynomial growth. 
Gromov proved the converse: polynomial growth implies virtually nilpotent. 

Gromov's paper is mind-blowing - the most astounding application of subRie- 
mannian geometry that I know of made by a human. (Cats and micro-organisms 
have made their own astounding applications.) Gromov scales the edges of the Cay- 
ley graph by e, then takes the limit as e — > 0. He proves, in essence, that the result 
converges to a Carnot group - a metric of subRiemannian type on a nilpotent Lie 
group - and from this the theorem easily follows. (I am stretching the truth here, 
but that is the spirit of Gromov's paper. There are many technicalities.) What 
I find so compelling about Gromov's paper is the going back and forth between 
the wonderful world of smooth metric spaces - Lie groups even - which I know 
and love, and the chopped up world of discrete objects that I find so frightening at 
times. Can we similarly go back and forth in dynamics? That is what I would like 
to see in some 'Kepler problem on a lattice.' 
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